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PREFACE 


During  the  last  three  years  considerable  progress  has  been  made  in  the  prediction  of 
unsteady  airloads  on  profiles  in  transonic  flow,  but  prediction  is  still  limited  to  potential 
inviscid  flow.  New  analytical  developments,  and  comparison  with  experimental  results,  need 
the  introduction  of  the  coupling  between  the  inviscid  How  and  the  boundary  layer.  Work  in 
this  field  is  at  its  beginning,  but  seems  very  promising  and  is  being  approached  in  many 
research  laboratories.  For  these  reasons,  the  Sub-Committee  on  Aeroelasticity  of  the 
Structures  and  Materials  Panel  considered  it  to  be  an  appropriate  time  for  the  presentation 
of  a pilot  paper  to  assist  in  the  preparation  of  a Specialists’  Meeting  in  the  Fall  of  1980. 

Dr  Yates,  in  this  publication,  proposes  an  original  approach  to  the  problem  which 
permits  a first  calculation  of  boundary  layer  effects.  He  derives  a new  kernel  of  the  lifting 
surface  integral  equation  which  takes  into  account  shear  flow  effects  in  the  boundary  layer. 
A very  lively  discussion  took  place  after  the  presentation  which  was  appreciated  as  an 
important  step  in  the  preparation  of  new  Panel  activities. 


G.COUPRY 

Chairman,  Sub-Committee 
on  Aeroelasticity 


v 


79  04  03  026 


CONTENTS 


Page 

SUMMARY  1 

LIST  OF  SYMBOLS  1 

I.  INTRODUCTION  1 

II.  CONCEPT  OF  VISCOUS  THIN  AIRFOIL  THEORY  3 

III.  THE  FLAT  PLATE  AIRFOIL  3 

IV.  DERIVATION  OF  AN  EQUIVALENT  INTEGRAL  EQUATION  4 

V.  NUMERICAL  RESULTS  - 6 

VI.  CONCLUSIONS  9 

FIGURES  10 


iv 


UNSTEADY  VISCOUS  THIN  AIRFOIL  THEORY 


John  E.  Yates 
Senior  Consultant 

Aeronautical  Research  Associates  of  Princeton,  Inc. 
50  Washington  Road,  P.0.  Box  2229 
Princeton,  New  Jersey  08540 


SUMMARY 

The  concept  of  viscous  thin  airfoil  theory  introduced 
in  Ref.  1 is  formulated  for  unsteady  incompressible 
flow.  The  theory  it'  developed  in  detail  for  a flat 
plate  airfoil  with  no  thickness  boundary  layer.  It  is 
shown  that  the  viscous  pressure-downwash  kernel  function 
has  a logarithmic  singularity  in  contrast  to  the  Cauchy 
singularity  of  inviscid  theory.  There  is  no  eigensolu- 
tion  and  no  need  for  a Kutta  condition  to  obtain  a unique 
solution.  It  is  shown  by  direct  numerical  solution  that 
for  Reynolds  number  greater  than  1000  the  viscous  and 
inviscid  results  are  virtually  the  same  except  in  the 
immediate  vicinity  of  the  trailing  edge.  There,  the 
pressure  loading  is  greater  than  inviscid  theory  would 
indicate  and  the  phase  of  the  complex  loading  is  less 
than  inviscid  theory.  The  results  are  in  qualitative 
agreement  with  experimental  results  of  Satyanarayana  and 
Davis,  Ref.  2.  The  effect  of  edge  bluntness  is  illu- 
strated for  the  case  of  steady  flow. 
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LIST  OF  SYMBOLS 

see  Eqs.  (5.1)  and  (5.2) 

airfoil  chord 

see  Eqs.  (5.2)  and  (5.3) 

unperturbed  and  perturbed  airfoil  surface  see  Fig.  1 

function  used  to  define  the  normal  displacement  of  the  airfoil  surface 

normalized  surface  displacement 

see  Eq.  (5.4) 

Bessel  function 

2^-  reduced  frequency  (based  on  semi-chord) 

modified  Bessel  function 

kernel  function,  see  Eq.  (4.19) 

integrated  kernel  function,  see  Eq.  (4.20) 

lift  distribution,  see  Eq.  (4.8) 

unit  vector  normal  to  airfoil  surface 

dimensionless  perturbation  pressure 

conjugate  pressure 

pressure  of  unperturbed  airfoil 

freestream  velocity 
v„c 

Reynolds  number 

time 

unit  vector  tangent  to  airfoil  surface 
Chebyshev  polynomial  of  the  first  kind 
Chebyshev  polynomial  of  the  second  kind 
potential  flowfield  of  unperturbed  airfoil 


v®  unperturbed  surface  speed 

V dimensionless  perturbation  velocity 

x , y Cartesian  coordinates,  see  Fig.  1 

y Eulers  constant  (0.57722) 

v kinematic  viscosity 

freestream  density 
a Re/4 

•(*  * ¥)* 

ill  stream  function 


Special  notation 

center  of  pressure 

lift  coefficient 

lift  curve  slope  (per  degree) 

moment  coefficient  (mid-chord  reference) 

2 

grad,  V vector  notation  for  divergence,  gradient  and  Laplacian  operators 

3t  + vo  * 8rad 
denotes  perturbation  quantity 
Im  real  and  imaginary  parts  of  a complex  quantity 


I . INTRODUCTION 

The  aeroelastlc  behavior  of  subsonic  aircraft  wings  can  be  calculated  with  well 
established  state  of  the  art  aerodynamic  methods.  Flutter  speeds  for  practical  aircraft 
wings  are  such  that  the  reduced  frequency  is  the  order  of  a few  tenths  and  the  conven- 
tional unsteady  aerodynamic  analyses  are  adequate.  We  remark,  however,  that  it  is  de- 
sirable to  modify  the  steady  state  aerodynamic  derivatives  to  agree  with  the  "best 
available"  steady  loads  data  (Ref.  3).  On  the  other  hand,  in  hydroelastic  applications, 
the  reduced  frequency  can  be  of  order  unityand  conventional  unsteady  aerodynamic  methods 
are  not  always  adequate.  For  example,  flutter  speed  predictions  for  hydrofoils  are  at 
best  unreliable  (Refs.  4 and  5). 

Several  authors  have  attempted  to  explain  some  of  the  discrepancies  between  theory 
and  experiment  (Ref.  6) . It  is  not  our  intent  to  give  a comprehensive  review  of  the 
aeroelastlc  work  but  to  call  attention  to  some  of  the  suspected  deficiencies  in  the 
aerodynamic  theory.  A significant  improvement  in  flutter  speed  prediction  was  obtained 
by  Carson  Yates  (Ref.  3)  when  steady  loads  data  were  used  to  correct  the  steady  state 
aerodynamic  derivatives.  These  results  suggest  a three-dimensional  tip  effect  that  is 
not  accounted  for  in  the  standard  linearized  aerodynamic  theory. 

Other  authors  (Refs.  7 and  8)  have  focused  on  the  Kutta  condition,  or  rather  its 
failure  at  higher  reduced  frequencies,  as  a prime  suspect  in  the  hydrofoil  problem. 
Various  suggestions  have  been  made  for  emperically  modifying  the  Kutta  condition  at 
higher  reduced  frequency.  Also,  several  authors  have  inquired  into  the  physical  basis 
of  the  Kutta  condition  by  including  viscosity  in  the  aerodynamic  theory.  This  was  ac- 
complished by  adopting  the  low  Reynolds  number  Incompressible  model  of  Oseen  to  calculate 
the  unsteady  aerodynamic  loads  on  a flat  plate  airfoil.  The  primary  conclusion  of  Shen 
and  Crimi  (Ref.  8)  was  that  the  viscosity  uniquely  determines  the  circulatory  part  of 
the  unsteady  load.  Furthermore,  the  loading  is  virtually  the  same  as  inviscid  theory 
with  Kutta  condition  with  small  corrections  of  order  (1//Ke)  . 

It  is  the  purpose  of  this  paper  to  re-examine  and  re-interpret  the  unsteady  viscous 
problem  in  the  light  of  more  recent  measurements  of  the  unsteady  pressure  distributions 
on  oscillating  airfoils,  and  also  to  compare  with  some  older  steady  state  lift  data. 

Some  preliminary  work  on  the  steady  state  problem  was  reported  in  Ref.  1.  Herein,  we 
formulate  the  problem  in  terms  of  an  integral  equation  (incompressible  Possio)  for  the 
ressure.  The  viscous  kernel  is  discussed  in  detail.  Numerical  results  are  presented 
or  the  flat  plate  airfoil  at  angle  of  attack,  oscillating  as  a rigid  body  and  subjected 
to  a sinusoidal  gust.  Our  results  corroborate  the  general  conclusions  of  previous 
authors  who  used  the  Oseen  model  but  are  re-interpreted  for  high  Reynolds  number  turbu- 
lent flow. 
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II.  CONCEPT  OF  VISCOUS  THIN  AIRFOIL  THEORY 

Consider  a Chin  two-dimensional  symmetric  airfoil  in  an  incompressible  viscous  flow 
at  zero  angle  of  attack  Fig.  la.  We  assume  that  the  Reynolds  number  is  very  large 
(>  10‘)  so  that  the  boundary  layer  is  turbulent  and  very  thin.  The  pressure  distribu- 
tion can  be  calculated  quite  accurately  with  potential  flow  for  the  basic  airfoil. 
Furthermore,  if  we  add  the  boundary  layer  displacement  thickness  to  the  geometric  thick- 
ness (dashed  line  in  Fig.  la)  we  can  pose  an  equivalent  potential  flow  problem.  We 
solve  for  the  potential  flow  on  the  mean  boundary  layer  corrected  airfoil  shape.  The 
surface  of^our^corrected  airfoil  (denoted  by  CQ  ) actually  moves  at  the  potential  flow 
velocity  t • vg  . The  mathematical  problem  can  be  posed  as  follows: 


dtv  vQ 
grad  vQ  + grad  PQ 
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(2.1) 


The  surface  speed  vg  is  to  be  calculated  from  the  solution  of  our  equivalent  potential 
flow  problem. 

Next,  we  suppose  that  the  airfoil  (now  understood  as  the  equivalent  airfoil)  surface 
is  perturbed  in  some  steady  or  unsteady  way  to  shape  C (Fig.  lb)  that  in  general  will 
produce  lift.  The  perturbation  flow  (assumed  to  be  small)  satisfies  the  equations 


where  f is  the  local  normal  displacement  of  the  surface. 

The  problem  of  "viscous  thin  airfoil  theory"  is  to  solve  the  viscous  perturbation, 

Eqs . (2.2)  subject  to  the  boundary  conditions,  Eq.  (2.5)  and  Eq.  (2.6),  and  the  infinity 
condition,  Eq.  (2.4).  If  we  drop  the  viscous  term  in  Eq.  (2.2)  and  the  perturbation  no 
slip  boundary  condition,  Eq.  (2.5),  the  .usual  inviscid  thin  airfoil  problem  is  recovered. 
However,  some  auxiliary  condition  (e.g.,  Kutta)  must  be  used  to  replace  the  missing  vis- 
cous terms.  Finally,  we  remark  that  an  interesting  and  practical  variation  of  the  fore- 
going development  is  to  replace  the  viscosity  by  the  local  eddy  viscosity  in  Eq.  (2.2). 

III.  THE  FLAT  PLATE  AIRFOIL 

For  the  remainder  of  this  report  we  restrict  our  attention  to  the  flat  plate  airfoil. 
Furthermore,  we  neglect  the  mean  boundary  layer  displacement  correction  in  which  case 


tv. 


(3.1) 


We  choose  the  semi-chord  c/2  as  the  unit  of  length,  the  frees tream  v«,  as  the  unit  of 
velocity  and  p„v&  as  the  unit  of  pressure.  The  dimensionless  viscous  thin  airfoil 
problem  for  the  flat  plate  becomes 

div  \f’  - 0 

Dp  J , 


(3.2) 

U(x,0±,t)  - 0 

(3.3) 

V(x,0±,t)  - -1  < x < 1 

(3.4) 

0 at  <» 

(3.5) 
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The  boundary 

condition , 

Eq.  (2.3), 

can  be  transferred  to  the  mean  surface  C 
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where 


D 3.3 
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(3.6) 

(3.7) 


and  F(x,t)  la  the  deflection  mode  of  the  airfoil.  Real  quantities  may  be  calculated 
with  the  relations, 


Time 

Distance 

Velocity 

Pressure 


(3.8) 
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The  momentum  equation,  Eq.  (3.2),  is  formally  identical  to  the  Oseen  equation. 
However,  the  justification  for  the  use  of  the  Eq.  (3.2)  for  high  Reynolds  number  is  com- 
pletely different  from  the  low  Reynolds  number  arguments  that  lead  to  the  Oseen  model. 

We  have  removed  the  symmetric  boundary  layer.  Also,  the  perturbation  viscous  boundary 
condition,  Eq.  (3.3),  is  different  from  the  full  no-slip  boundary  condition  used  with 
the  Oseen  equation.  The  basic  theory  differs  even  further  from  the  Oseen  model  if  we 
use  a Reynolds  number,  Eq.  (3.7),  based  on  a typical  eddy  viscosity.  Many  of  the  results 
we  derive  have  a formal  similarity  to  the  results  of  Chu  (Ref.  7)  and  Shen  and  Crimi 
(Ref.  8)  who  used  Oseen's  equation.  However,  the  interpretation  of  our  results  is  for 
high  Reynolds  number.  Also,  the  theory  can  be  extended  in  a logical  self  consistent  way 
to  account  for  geometric  and  boundary  layer  displacement  thickness. 

IV.  DERIVATION  OF  AN  EQUIVALENT  INTEGRAL  EQUATION 

It  is  convenient  to  recast  the  viscous  thin  airfoil  problem  in  the  form  of  an  inte- 
gral equation.  The  procedure  is  much  the  same  as  with  inviscid  theory.  For  simplicity 
we  consider  simple  harmonic  motion  and  derive  a fundamental  solution  corresponding  to  a 
pressure  dipole  located  at  the  origin.  By  superposition  of  fundamental  solutions  we  can 
write  down  the  analog  of  Possio's  integral  equation  for  incompressible  viscous  flow. 

Q can  be  expressed  in  the  form 

(4.1) 


(4.2) 


is  the  usual  reduced  frequency  based  on  the  semi-chord.  We  retain  the  same  symbol  Q 
for  the  complex  amplitude  of  the  dependent  variables.  The  convective  operator,  Eq.  (3.6) 
becomes 


For  simple  harmonic  motion  any  dependent  variable 
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Next,  we  introduce  a streamfunction,  41  , such  that 


(4.3) 


$ = grad  p x it 


(4.4) 


satisfies  the  continuity  equation.  The  momentum  equation,  Eq.  (3.2),  can  be  expressed 
in  the  simple  form 


grad  P - grad  P * k 

where  P is  the  conjugate  pressure  (Ref.  7) 

* 


(4.5) 


(4.6) 


Both  P and  P are  solutions  of  Laplace's  equation  and  Eq.  (4.5)  is  the  vector  form 
of  the  Cauchy  Rieman  equations.  Thus, 


P = P + iP* 

is  an  analytic  function  of  the  complex  cariable  z = x + iy  . Below,  we  choose  to  work 
with  the  real  functions  P and  P*  . However,  it  is  worth  noting  that  a complex  repre- 
sentation for  the  pressure  is  possible  even  with  viscous  effects  included.  A complex 
analytic  potential  does  not  exist  for  the  viscous  problem. 
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By  symmetry  of  the  boundary  conditions  the  pressure  can  be  expressed  as  a super- 
position of  dipoles  along  the  airfoil  cho^rd;  i.e.. 


in  R dC 


L(x)  - P(x,0”)  - P(x,0+) 

R - /(x  - O2  + yZ 

One  can  show  by  direct  substitution  into  Eq.  (4.5)  that 

.1 


in  R dC 


(4.10) 


The  streamfunction  corresponding  to  our  assumed  pressure  distribution  is  to  be 
found  by  solution  of  Eq.  (4.6).  Consider  the  equation 


27  v * ~ B*  37  in  R 

In  the  absence  of  viscosity  the  streamfunction  is  given  by 


(4.11) 


37  ln  R dc 


so  that 


vV°> 


2"  4 [e 


H(x) 6 (y) 


(4.12) 


(4.13) 


where  H(x)  and  6(y)  are  the  Heaviside  step  and  Dirac  delta  functions  respectively. 

The  inviscid  solution  is  a periodic  vortex  sheet  that  is  shed  from  the  dipole  and  extends 
unattenuated  downstream  to  infinity.  The  effect  of  viscosity  is  to  smooth  out  the  dis- 
continuity and  cause  the  wake  to  decay. 


Next,  we  write  the  solution  of  Eq.  (4.11)  in  the  form 

* - *<°>  + *<!> 


where 


is  vV1} 


d«p _ it  3 r _ikx  w/  n 
= H(x)6(y)J 


(4.14) 


(4.15) 


The  solution  can  easily  be  expressed  in  a form  similar  to  Eq.  (4.12).  The  final  form  of 
ip  is 


a(x- 

in  R + e 


Kq(o  R)  d£ 


(4.16) 


(4.17) 


and  Kq  is  the  modified  Bessel  function. 

The  perturbation  velocity  field  corresponding  to  Eq.  (4.15)  is  everywhere  continuous. 
It  follows  by  antisymnetry  that  the  x-component  of  velocity  is  zero  on  the  x-axis.  Thus, 
the  no-slip  boundary  condition,  Eq.  (3.3),  is  satisfied  for  an  arbitrary  line  distribu- 
tion of  fundamental  solutions.  Also,  the  perturbation  velocity  decays  at  infinity  so 
that  Eq.  (3.5)  is  satisfied.  The  final  step  is  to  superimpose  the  fundamental  solutions 
and  apply  the  boundary  condition,  Eq.  (3.4),  on  the  normal  component  of  velocity.  The 
final  form  of  the  integral  equation  for  L(x)  becomes 


with  the  kernel  defined  by 
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We  refer  to  K(x)  as  the  kernel  and  K^(x)  as  the  integrated  kernel.  The  mathematical 


properties  of  K(x)  determine  the  basic  character  of  the  solution  of  Eq . (4.18)  and  are 
discussed  below.  For  numerical  purposes  it  is  more  convenient  to  evaluate  the  inte- 
grated kernel. 


Discussion  of  the  Kernel 


The  interpretation  of  the  kernel  is  the  same  in  viscous  thin  airfoil  theory  as  it 
is  in  classical  inviscid  theory;  i.e.,  it  is  the  downwash  velocity  induced  by  a dipole 
pressure  excitation  at  the  origin.  For  high  Reynolds  number  the  most  important  effect 
of  viscosity  is  near  the  origin.  We  can  see  this  easily  in  the  case  of  steady  flow 
(Ref.  1)  where 
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The  viscous  kernel  has  a weak  symmetric  logarithmic  singularity  while  the  inviscid  ker- 
nel has  a Cauchy  singularity.  This  statement  is  also  true  for  the  unsteady  kernel. 

From  the  theory  of  integral  equations  of  the  first  kind  (see  Ref.  9)  it  is  known  that 
Eq.  (4.18)  has  a unique  solution  for  a logarithmic  kernel  and  has  an  eigensolution  for 
a Cauchy  kernel.  In  inviscid  theory  we  fix  the  magnitude  of  the  eigensolution  by  invok- 
ing the  Kutta  condition  at  the  trailing  edge.  In  viscous  thin  airfoil  theory  no  auxil- 
iary condition  is  required  to  solve  the  integral  equation. 


Another  important  property  of  the  kernel  is  the  difference  between  upstream  and 
downstream  influence.  In  steady  flow,  for  example,  the  inviscid  kernel  is  symmetric 
upstream  and  downstream.  Viscosity  causes  a reduction  in  the  downwash  downstream  of  the 
load.  Analytically  the  kernel  is  approximately 
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(4.22) 


The  steady  state  kernel  is  plotted  in  Fig.  2 as  a function  of  ox  . The  nature  of  the 
singularity  and  the  enhanced  downstream  influence  are  clearly  evident  in  the  viscous 
kernel . 


The  unsteady  viscous  kernel  exhibits  another  fundamental  difference  from  its  invis- 
cid counterpart.  The  inviscid  vortex  wake  (downwash  field)  does  not  decay  while  the 
viscous  wake  decays  exponentially  like 


-k2x/2o 


(4.23) 


(see  Fig.  3).  This  feature  of  the  kernel  is  physically  appealing  in  that  a real  wake 
does  decay  ultimately  due  to  turbulence.  If  we  base  the  Reynolds  number  on  a typical 
wake  eddy  viscosity  then  the  decay  formula,  Eq . (4.23)  becomes  even  more  meaningful. 

The  effect  of  Reynolds  number  on  the  integrated  unsteady  kernel  near  the  origin  is  illu- 
strated in  Fig.  4.  For  fixed  x > 0 the  difference  between  the  viscous  and  inviscid 
kernel  vanishes  like  1//Ee  . Upstream,  the  convergence  is  exponential. 


V.  NUMERICAL  RESULTS 


To  solve  the  basic  integral  equation  (4.18),  we  expand  the  load  distribution  L(x) 
in  a Chebyshev  series;  i.e., 

N 

X T (x) 

L(x)  = > A„  -^-2 , (5.1) 


where  Tn(x)  is  the  Chebyshev  polynomial  of  the  first  kind  (Ref.  10).  Substitute  Eq. 
(5.1)  into  Eq.  (4.18),  multiply  by  Um(x)  (Chebyshev  polynomial  of  the  second  kind)  and 
integrate  form  -1  to  1 . The  result  is  a matrix  equation  for  the  unknown  coefficients, 
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Three  special  cases  of  the  mode  shape,  downwash  and  the  coefficients  e are  noted 
below: 
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where  Jm  is  the  Bessel  function  of  the  first  kind.  The  coefficients  are  functions 

of  the  reduced  frequency  and  the  Reynolds  number  and  must  be  evaluated  numerically.  We 
remark  that  it  is  convenient  to  introduce  the  trigonometric  substitution  in  Eq.  (5.3) 


x “ cos  6 


5 = cos 


(5.8) 


in  which  case 
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(5.9) 


The  integrated  viscous  kernel  is  regular  at  the  origin  so  that  the  integrals  in  Eq. 

(5.9)  can  be  evaluated  by  any  standard  quadrature  routine.  While  the  general  theory  of 
integral  equations  of  the  airfoil  type  is  not  our  main  concern  in  this  paper,  we  point 
out  that  Williams  (Ref.  11)  has  derived  a general  solution  of  Eq.  (4.18)  in  terms  of  two 
fundamental  solutions.  The  latter  correspond  essentially  to  the  cases  of  translation 
and  pitch  about  mid-chord.  Also,  Williams  points  out  (Ref.  12)  that  the  evaluation  of 
matrix  coefficients  like  Eq.  (5.9)  can  be  reduced  to  a single  quadrature  over  the  Fourier 
transform  of  the  kernel  function.  This  might  be  a more  practical  method  than  the  direct 
evaluation  of  the  kernel. 


For  any  finite  N and  in  the  absence  of  viscosity,  the  matrix  C^n  is  singular  and 
has  the  rank  N-l  . There  is  one  eigensolution  whose  amplitude  is  fixed  by  invoking  the 
Kutta  condition  at  the  trailing  edge.  In  practice  the  Kutta  condition  is  satisfied  by 
choosing  a set  of  pressure  mode  shapes  that  individually  vanish  at  the  trailing  edge. 

In  Eq.  (5.1)  we  have  allowed  for  the  possibility  of  a square  root  singularity  in  the  pres- 
sure loading  at  the  leading  and  trailing  edges.  With  the  viscous  kernel  the  matrix  Cjm 
is  not  singular  and  there  is  no  eigensolution.  The  circulatory  component  of  the  unsteady 
loading  is  uniquely  determined  for  any  finite  Reynolds  number.  We  are  of  course  rephrasing 
in  matrix  terms  our  earlier  observation  that  the  integral  equation,  Eq.  (4.18),  has  a 
unique  solution  for  the  viscous  kernel  because  of  the  logarithmic  singularity. 

It  is  a curious  theoretical  fact  (Ref.  9)  that  the  inclusion  of  viscosity  in  the 
analysis  does  not  alter  the  square  root  edge  singularity  in  the  pressure.  Only  by  in- 
cluding the  nonlinear  effect  of  edge  bluntness  can  this  singularity  be  removed. 


The 
terms  of 


lift  and  moment  coefficients  and  center  of  pressure  are  readily  calculated  in 
the  Chebyshev  coefficients;  i.e., 
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In  Fig.  5,  the  flat  plate  lift  curve  slope  calculated  from  Eq.  (4.18)  is  plotted  as 
a function  of  Reynolds  number  and  compared  with  the  limit  case  of  inviscid  theory  plus 
the  Kutta  condition.  The  viscous  flat  plate  solution  converges  to  the  inviscid  limit  as 
1//Ke  . For  Re  > 1000  , the  viscous  and  inviscid  results  are  virtually  the  same.  These 
numerical  results  are  in  agreement  with  the  conclusion  of  Shen  and  Criml  (Ref.  8).  Even 
though  we  have  allowed  for  the  possibility  of  a square  root  singularity  in  the  loading 
at  the  trailing  edge,  the  effect  of  viscosity  is  to  remove  the  singularity  in  accordance 
with  the  empirical  Kutta  condition. 

The  insight  that  the  viscous  analysis  provides  is  that  the  "Kutta  condition"  is  a 
consequence  of  the  integrated  effect  of  viscosity  over  the  entire  wetted  surface  of  the 
airfoil.  This  point  can  be  made  stronger  in  the  case  of  unsteady  flow.  For  example, 
in  Figs.  6a  and  6b  we  compare  our  theory  with  the  experimental  results  of  Satyanarayana 
and  Davis  (Ref.  2) . There  is  no  difference  between  the  viscous  and  inviscid  theory 
except  over  roughly  the  last  5%  of  the  airfoil.  There  the  magnitude  of  the  load  is  greater 
than  the  inviscid  result  that  drops  to  zero  as  the  square  root  of  the  distance  from  the 
trailing  edge  in  accordance  with  the  Kutta  condition.  Also,  there  is  a drop  in  the  phase 
of  the  loading.  These  results  are  in  qualitative  agreement  with  the  experimental  results. 
However,  the  experiments  indicate  much  greater  deviations  from  local  Kutta  behavior  in 
particular  at  the  reduced  frequency  of  1.23  (see  Fig.  6b). 

The  main  point  is  that  the  loading  over  the  majority  of  the  frontal  portion  of  the 
airfoil  can  be  calculated  quite  accurately  with  inviscid  theory  plus  the  Kutta  condition 
even  though  the  local  behavior  of  the  loading  at  the  trailing  edge  may  in  no  way  resemble 
the  potential  flow.  It  Is  the  integrated  effect  of  viscosity  that  is  responsible  for  the 
circulatory  component  of  the  steady  or  unsteady  loading.  This  is  the  reason  why  local 
modifications  (e.g.,  bluntness)  to  an  airfoil  trailing  edge  can  be  made  without  making 
gross  changes  in  its  overall  lift  performance  (Ref.  13). 

A final  result  to  illustrate  our  main  point  is  shown  in  Fig.  7 where  we  subject  an 
airfoil  to  a sinusoidal  gust  at  a reduced  frequency  of  2.0.  Again,  the  viscous  and 
inviscid  loadings  are  indistinguishable  over  the  front  portion  of  the  airfoil.  Near  the 
trailing  edge  the  load  increases  and  the  phase  dips  below  its  inviscid  counterpart. 

Recent  measurements  of  unsteady  pressure  due  to  buffeting  (Ref.  14)  indicate  that  the 
enhanced  trailing  edge  loading  does  occur. 

Effect  of  Edge  Bluntness 

In  Ref.  1,  we  presented  some  results  for  steady  flow  that  illustrate  the  importance 
of  leading  edge  bluntness.  The  measured  lift  curve  slope  of  a large  collection  of  stan- 
dard airfoils  (Ref.  15)  is  15  to  207.  less  than  the  value  predicted  by  inviscid  theory 
with  Kutta  condition.  Typical  experimental  results  for  Reynolds  numbers  between  10“  and 
10'  are  shown  in  Fig.  8.  The  theoretical  result  was  calculated  with  an  ad  hoc  modifica- 
tion of  viscous  thin  airfoil  theory  that  in  essence  suppresses  the  effect  of  viscosity 
in  the  immediate  vicinity  of  the  leading  edge.  The  argument  is  that  if  the  leading  edge 
were  blunt  like  an  actual  subsonic  airfoil,  then  the  square  root  leading  edge  singularity 
of  the  pressure  would  not  be  present.  The  effect  of  viscosity  combined  with  a sharp  lead- 
ing edge  is  to  grossly  overestimate  the  leading  edge  load  and  total  circulation.  Our 
ad  hoc  correction  suppresses  chis  singular  effect  and  yields  a formula  for  the  lift  curve 
slope  that  has  a logarithmic  Reynolds  number  correction.  The  order  of  magnitude  and 
trend  of  the  correction  is  in  very  good  agreement  with  experimental  data  as  shown  in  Fig. 

8.  Because  of  the  success  of  our  ad  hoc  correction,  considerable  emphasis  is  being  placed 
cn  a more  rational  correction  for  edge  bluntless  in  our  current  research  work. 

Comments  on  Asymptotic  Trailing  Edge  Theories 

In  recent  years  several  papers  (Refs.  16,  17,  18)  have  considered  the  detailed  flov 
around  an  airfoil  with  a mathematically  sharp  trailing  edge.  The  problem  is  to  match 
every  detail  an  inner  viscous  boundary  layer  flow  to  an  external  potential  flow.  Multi- 
layered asymptotic  analysis  is  necessary  to  complete  the  matching.  The  most  practical 
application  of  this  type  of  analysis  is  the  recent  strong  interaction  theory  of  Melnik 
(Ref.  18).  The  lift  and  drag  of  an  airfoil  in  steady  compressible  flow  with  an  attached 
turbulent  boundary  layer  can  now  be  calculated  with  great  accuracy. 

All  of  these  analyses  rely  heavily  on  a highly  detailed  asymptotic  analysis  of  the 
trailing  edge  region  to  connect  the  boundary  layers  on  the  upper  and  lower  surfaces. 

The  trailing  edge  must  be  mathematically  sharp  and  in  the  innermost  solution  near  the 
trailing  edge  the  "least  singular"  branch  of  the  possible  solutions  is  .-hosen.  This  is 
equivalent  to  applying  a "Kutta"  condition.  It  seems  ironical  that  an  empirical  mathe- 
matical statement  must  be  made  to  establish  the  uniqueness  of  such  a beautiful  asymtotic 
theory.  Yet,  we  can  illustrate  very  simply  why  such  a condition  must  b§  used  with  con- 
ventional boundary  layer  analysis.  If  we  replace  the  viscous  term  (\>VJv)  in  our  Eq.(2.2) 

with  its  "parabolized"  boundary  layer  counterpart^  -3— ) then  it  i3  not  difficult  to 

' 3y*  ' 

show  (we  omit  the  proof  here)  that  the  resulting  kernel  function  in  our  viscous  thin 
airfoil  theory  still  has  a Cauchy  singularity.  The  solution  for  the  lift  distribution 
will  have  an  eigensolution  and  some  auxiliary  criteria  (Kutta)  must  be  used  to  establish 
a unique  solution.  The  full  elliptic  behavior  of  the  viscous  term  must  be  retained  in 
the  Navier  Stokes  equations  if  the  need  for  a Kutta  condition  is  to  be  eliminated.  There 
is  no  reason  to  suspect  that  the  nonlinear  terms  in  the  fluid  equations  can  change  this 
requirement . 
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Finally,  we  remark  that  the  concept  and  spirit  of  viscous  thin  airfoil  theory  is 
completely  different  from  the  asymptotic  theories.  Our  basic  working  principle  is  that 
the  most  important  viscous  corrections  in  the  gross  aerodynamic  properties  needed  for 
aeroelastic  calculations  can  be  calculated  without  recourse  to  highly  detailed  boundary 
layer  analysis.  To  the  extent  that  the  Kutta  condition  is  in  itself  the  most  important 
viscous  correction,  the  principle  has  already  served  us  well. 

VI.  CONCLUSIONS 

We  have  formulated  the  concept  of  unsteady  viscous  thin  airfoil  theory  and  applied 
the  theory  to  a flat  plate  airfoil.  For  this  special  case  our  results  are  formally 
equivalent  to  those  of  Chu  (Ref.  7),  and  Shen  and  Crimi  (Ref.  8),  who  use  the  Oseen 
equation.  The  present  theory  is  valid  for  high  Reynolds  number  turbulent  flow  in  contrast 
to  the  Oseen  model. 

The  main  conclusion  of  our  work  is  that  viscosity  is  the  most  important  missing 
element  of  physics  in  our  potential  flow  theories.  Without  viscosity  one  must  invoke  an 
auxiliary  condition  such  as  the  Kutta  condition  to  determine  the  circulatory  component 
of  the  loading.  The  viscous  theory  indicates  that  the  development  of  circulation  is 
due  to  the  integrated  effect  of  viscosity  over  the  entire  wetted  surface  of  the  airfoil. 
Local  deviations  from  smooth  Kutta  flow  at  the  trailing  edge  are  quite  likely  to  occur 
on  any  real  airfoil  without  altering  the  overall  aerodynamic  loads  required  in  an  aero- 
elastic analysis. 


An  ad  hoc  correction  for  leading  edge  bluntness  is  presented  for  the  case  of  steady 
flow.  The  correction  to  the  lift  curve  slope  is  of  order  l/tn  Re  and  agrees  in  magnitude 
and  trend  with  high  Reynolds  number  steady  state  lift  curve  data.  We  conclude  that  lead- 
ing edge  bluntness  must  be  accounted  for  in  viscous  thin  airfoil  theory  to  avoid  a severe 
overestimation  of  the  leading  edge  load  and  total  circulation. 
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Figure  4.  Effect  of  Reynolds  number  on  the  integrated  kernel  near 


Figure  5.  Solution  of  the  viscous  airfoil  equation  for  an  airfoil  at  constant 
angle  of  attack. 
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Figure  6a.  Lift  distribution  (magnitude  and  phase)  on  an  airfoil  oscillating  about 
quarter  chord;  Re  = 1000,  k - 1.02  . 
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Figure  6b.  Lift  distribution  (magnitude  and  phase)  on  an  airfoil  oscillating  about 
quarter  chord;  Re  = 1000,  k = 1.23  . 
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Figure  7.  Lift  distribution  on  an  airfoil  in  a sinusoidal  gust;  Re  = 1000, 
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^.Abstract 

1 During  the  last  three  years  considerable  progress  has  been  made  in  the  prediction  of  unsteady 
airloads  on  profiles  in  transonic  flow,  but  prediction  is  still  limited  to  potential  in  viscid  flow. 
New  analytical  developments,  and  comparison  with  experimental  results,  need  the  introduc- 
tion of  the  coupling  between  the  inviscid  flow  and  the  boundary  layer.  Work  in  this  field  is 
at  its  beginning,  but  seems  very  promising  and  is  being  approached  in  many  research 
laboratories.  For  these  reasons,  the  Sub-Committee  on  Aeroelasticity  of  the  Structures  and 
Materials  Panel  considered  it  to  be  an  appropriate  time  for  the  presentation  of  a pilot  paper 
to  assist  in  the  preparation  of  a Specialists’  Meeting  in  the  Fall  of  1980. 

^ Dr  Yates,  in  this  publication,  proposes  an  original  approach  to  the  problem  which  permits  a 
first  calculation  of  boundary  layer  effects.  He  derives  a new  kernel  of  the  lifting  surface 
integral  equation  which  takes  into  account  shear  flow  effects  in  the  boundary  layer.^A  very 
lively  discussion  took  place  after  the  presentation  which  was  appreciated  as  an  important 
step  in  the  preparation  of  new  Panel  activities. 
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